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ABSTRACT 

The iron mass in galaxy clusters is about 6 times larger than could have been 
produced by core-collapse supernovae (SNe), assuming the stars in the cluster formed 
with a standard initial mass function (IMF). SNe la have been proposed as the alterna- 
tive dominant iron source. Different SN la progenitor models predict different "delay 
functions" , between the formation of a stellar population and the explosion of some of 
its members as SNe la. We use our previous measurements of the cluster SN la rate 
at high redshift to constrain SN la progenitor models and the star-formation epoch in 
clusters. The low observed rate of cluster SNe la at z ~ — 1 means that, if SNe la 
produced the observed amount of iron, they must have exploded at even higher z. 
This puts a > 95% upper limit on the mean SN la delay time of r < 2 Gyr (< 5 Gyr) 
if the stars in clusters formed at Zf < 2 (zf < 3), assuming H Q — 70 km s^ 1 Mpc -1 . 
In a companion paper, we show that, for some current versions of cosmic (field) star 
formation history (SFH), observations of field SNe la place a lower bound on the delay 
time, t > 3 Gyr. If these SFHs are confirmed, the entire range of r will be ruled out. 
Cluster enrichment by core-collapse SNe from a top-heavy IMF will then remain the 
only viable option. 
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1 INTRODUCTION 

Galaxy clusters, by virtue of their deep gravitational po- 
tentials, are "closed boxes", from which little matter can 
escape. They thus constitute ideal sites to study the time- 
integrated enrichment of the intergalactic medium (see e.g., 
Buote 2002, for a recent review). Observations of galaxy 
clusters have revealed a number of intriguing puzzles. One 
such puzzle follows from X-ray spectroscopy, and shows that 
the intracluster medium (ICM) gas consistently has a sur- 
prisingly high iron abundance, with a "canonical" value of 
about 0.3 the Solar abundance (e.g., Mushotzky & Loewen- 
stein 1997; Fukazawa et al. 1998; Finoguenov et al. 2000; 
White 2000). Combined with the fact that most of the bary- 
onic mass of the cluster is in the ICM, this translates into a 
large mass of iron. The tight correlation between total iron 
mass and stellar light from the early-type galaxy population 
(and the lack of a correlation with late-type galaxies) suggest 
that the stellar population, whose remnants and survivors 
now populate the early-type galaxies, produced the ICM 
iron (Renzini et al. 1993; Renzini 1997). The near constancy 
of the ICM iron abundance with cluster mass (Renzini 1997; 



Lin, Mohr, & Stanford 2003) and with redshift out to z ~ 1 
(Tozzi et al. 2003) further argues against a significant role 
in ICM iron enrichment for recent infall and disruption of 
metal-rich dwarf galaxies. However, the iron mass is at least 
several times larger than that expected from core-collapse 
supernovae (SNe), based on the present-day stellar masses, 
and assuming a standard stellar initial mass function (IMF; 
e.g., Renzini et al. 1993; Loewenstein 2000). The problem is 
aggravated if one considers the large mass of iron that exists 
in the cluster galaxies themselves. 

A possibly related problem is the energy budget of the 
ICM gas and the "entropy floor" observed in clusters, which 
suggest a non-gravitational energy source to the ICM, again 
several times larger than the expected energy input from 
core-collapse SNe (e.g., Lloyd-Davies, Ponman, & Cannon 
2000; Tozzi & Norman 2001; Brighenti & Mathews 2001; 
Pipino et al. 2002). 

Proposed solutions to these problems have included an 
IMF skewed toward high-mass stars (so that a large number 
of iron-enriching core-collapse SNe are produced per present- 
day unit stellar luminosity), or a dominant role for SNe la 
in the ICM Fe enrichment. For example, Brighenti & Math- 
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ews (1998) calculate that the cluster SN la rate must be 
> 4.8/i 2 SNu today and > 9.6/i 2 SNu at z = 1 to explain 
production of most of the iron in the ICM with SNe la. 
[1 SNu= 1 SN century-^lO 10 !^©)- 1 .] This contrasts with 
the local elliptical-galaxy SN rate of (0.28 ± 0.12)h 2 SNu 
(Cappellaro et al. 1999) and argues against the type-la en- 
richment scenario. On the other hand, Renzini (1997) has 
pointed out that the approximately Solar abundance ratios 
measured in cluster ellipticals argue for a mix of SN-types, 
and hence also an IMF, that are not too different from the 
mix that exists in the Milky Way. Attempts to derive the 
SN mix in clusters by means of direct X-ray measurements 
of element abundance ratios in the ICM are still ambigu- 
ous, with some results favoring a dominance of la's (e.g., 
Buote 2002; Tamura et al. 2002) and others a dominance of 
core-collapse SNe (e.g., Lima Neto et al. 2003; Finoguenov, 
Burkert, & Boehringer 2003). 



In Gal- Yam, Maoz, & Sharon (2002), we used multi- 
ple deep Hubble Space Telescope (HST) archival images of 
galaxy clusters to discover distant field and cluster SNe. The 
sample was composed of rich clusters, with X-ray temper- 
atures in the range 4 — 12 keV, and a median of 9 keV. 
The candidate type-la SNe in the clusters then led to 
an estimate of the SN la rate in the central 250/i -1 kpc 
of medium-redshift (0.18 < z < 0.37) and high-redshift 
(0.83 < z < 1.27) cluster sub-samples. The measured rates 
are low. To within errors, they are not different from SN la 
measurements in field environments, both locally (Cappel- 
laro et al. 1999) and at high redshift (Pain et al. 2002; Tonry 
et al. 2003). It was argued that our 95 per cent upper lim- 
its on the cluster SN la rates rule out the particular model 
by Brighenti & Mathews (1998) for SNe la as the primary 
source of iron in the ICM. 



The issue of SN la rate vs cosmic time is closely tied 
to the presently unsolved question regarding the progeni- 
tor populations of SNe la. Different models predict different 
delay times between the formation of a stellar population 
and the explosion of some of its members as SNe la (e.g., 
Ruiz-Lapuente & Canal 1998; Yungelson & Livio 2000, and 
references therein). The SN la rate vs time in a given envi- 
ronment (e.g., cluster, or field) will then be a convolution of 
the star-formation history in that environment with a "de- 
lay" or "transfer function" , which is the SN la rate vs time 
following a brief burst of star formation [e.g., Sadat et al. 
1998; Madau, Delia Valle, & Panagia 1998 (MDP); Dahlen 
& Fransson 1999; Sullivan et al. 2000]. 



In the present paper, the cluster iron mass problem is re- 
visited, and some simple relations connecting the iron mass 
and the SN rate to updated values of the various observables 
are derived. We then use the observed upper limits on the 
cluster SN la rate to set constraints on the SN la progenitor 
models, on the formation of stellar populations in galaxy 
clusters, and on the cluster enrichment scenario. Through- 
out the paper we assume a flat cosmology with Q m = 0.3 
and Qa = 0.7, and a Hubble parameter of Ho = 70 km s _1 
Mpc" 1 . 



2 THE IRON PROBLEM, AND THE SN 
RATES NEEDED TO RESOLVE IT 

Detailed models of cluster metal enrichment have been cal- 
culated previously (e.g., Buote 2002; Pipino et al. 2002; 
Finoguenov et al. 2003). However, the specific problem of 
the total iron mass in clusters can be formulated rather sim- 
ply as a function of observational parameters. The ratio of 
observed iron mass to iron mass expected from core-collapse 
SNe is 



observed 

M /bar (fgas ^Fe-gas + /* Z^e*) 

Mfc-snii ~ M/ bar /» /(> 8 Mo) /fc-snii ' 



(1) 



Here M is the mass of a cluster, /bar is the baryon mass 
fraction, f gas is the mass fraction of the baryons in the ICM, 
and ZF C _gas is the mass fraction of that gas in iron. Similarly, 
/* = 1 — fgas is the mass fraction of the baryons in stars and 
Z Fe * is the iron abundance of the stars. In the denominator, 
/(> 8Mq) is the ratio of the initial stellar mass in stars of 
mass > 8Mq (i.e., stars that underwent core-collapse), to 
the mass in stars of lower mass, and / Fc -snii is the iron 
yield of core-collapse SNe, expressed as a fraction of the 
progenitor masses. 

Lin et al. (2003) have recently derived improved esti- 
mates of stellar and gas mass fractions in clusters using in- 
frared data from the 2MASS survey. For rich clusters they 
find f gas — 0.9 and /» = 0.1. Their stellar mass fraction is 
based on the 2MASS K-band luminosity function (Kochanek 
et al. 2001) combined with dynamical stellar mass-to-light 
ratio measurements by Gerhard et al. (2001). Ettori (2003) 
has recently suggested that a significant fraction, 6-38%, of 
cluster baryons may be in a yet-undetected warm gas. How- 
ever, if this new component has a similar iron abundance 
to that of the hot ICM, our arithmetic will not be affected. 
We adopt the "canonical" ICM iron abundance in rich clus- 
ters of Z Fo _ g as ~ 0.3 ZoFe (e.g., Mushotzky & Loewenstein 
1997; Fukazawa et al. 1998; Finoguenov et al. 2000; White 
2000). This iron abundance relates to a photospheric Solar 
iron mass abundance of Zq Fc = 0.0026 found by Anders 
& Grevesse (1989). [Using the updated photospheric Solar 
abundance of Zqfo — 0.00177, given by Grevesse & Sauval 
(1999), which also agrees with the meteoritic Solar value of 
Anders & Grevesse (1989), would imply simply raising the 
ICM value accordingly]. Most of the stellar mass is in the 
elliptical galaxies, for which we adopt Z Fc » ~ 1.2 Zqf s , the 
median found by j0rgensen (1999) for early- type galaxies in 
Coma. 

To estimate, /(> 8M ), the ratio of exploding to non- 
exploding initial stellar masses, an IMF, dN '/dm, must be 
assumed: 

f™ P dN/dm m dm 



1™*° dN/dm m dm 



(2) 



where mi ow and m up are the lower and upper mass cutoffs 
of the IMF. For a Salpeter (1955) IMF, dN/dm oc m~ 2 ' 35 . 
Baldry and Glazebrook (2003) have recently modeled the 
local UV-to-IR luminosity density of galaxies assuming a 
range of IMFs and SFHs, and found the data to be consistent 
with the Salpeter IMF. A Salpeter slope with mi ow — O.IMq 
and m up = 100M©, gives /(> 8M Q ) w 0.16. However, other 
IMFs have been proposed. Figure 1 shows the dependence 
of /(> 8M ) on m low for single-power-law IMFs (such as 
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Combining the estimates above, we can parametrize the 
iron problem as 
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Figure 1. The ratio of the initial stellar mass in stars of mass 
> 8Mq (i.e., stars that underwent core-collapse), to the initial 
mass in stars of lower mass (i.e., the present-day cluster stellar 
mass), shown as a function of the lower mass cutoff of the IMF, 
for various IMFs, as marked. The dotted line shows the value 
of /(> 8Mq) needed to produce the observed iron mass in rich 
clusters with core-collapse SNe only, for the fiducial values of 
the other parameters in Eq. 3. Single power-law IMFs (including 
Salpctcr 1955) are marked with their index. The "standard" non- 
top-heavy IMFs are from Kennicutt (1983; a = —1.4, —2.5 in 
the mass ranges 0.1 — IMq, 1 — IOOMq, respectively); Gould ct 
al. (1997; a = -0.9, -2.21, -2.35 for 0.1 - 0.6M Q , 0.6 - IMq, 
1 - 100M Q ); Scalo (1998; a = -1.5, -2.7, -2.3 for 0.1 - IMq, 
1 - 10M Q , 10 - 100M Q ); and Kroupa (2000; a = -1.3, -2.3 for 
0.1 - O.5M , 0.5 - IOOMq). 



Salpeter's) of various indices a, and for "standard" IMFs. 
Varying m up has a weak effect on /(> 8Mq), as long as the 
IMF is steep enough. 

Core collapse models generally agree on an iron yield 
of about O.IMq per SN (e.g., Thielemann, Nomoto, & 
Hashimoto 1996). In observed core-collapse SNe, estimates 
of the total yield of Ni 56 (which decays to Fe 56 and drives the 
optical luminosity of SNe) have been obtained by different 
methods - the luminosity of the radioactive tail, the lumi- 
nosity of the plateau phase, and the Ha luminosity in the 
nebular phase (Elmhamdi, Chugai, & Danziger 2003). The 
different methods give consistent results for a given SN, but 
there is a large scatter among different SNe (e.g., Zampieri 
et al. 2003) and a range of yields as large as 0.0016 to 0.26 
M Q (Hamuy 2003), with a mean of 0.05M© (Elmhamdi et 
al. 2003). Progenitor masses are more difficult to estimate, 
but there are some indications that the iron yield and the 
progenitor mass are correlated, i.e., the iron yield is some 
fraction of the progenitor mass, of order 0.5-1% (P. Maz- 
zali, private communication). Since the ratio of the mean 
iron yield found by Elmhamdi et al. (2003; 0.05M©) to the 
minimum progenitor mass (8Mq) is 0.063%, and the ratio 
is probably even smaller for those events with more massive 
progenitors, we will err conservatively by assuming a large 
fractional iron mass yield of /f c snii ~ 0.01. 



Fc — observed 



= 6.3 



/(> 8M ) 
0.16 



{fgas Z~Fe — gas ~\~ f* ^Fe*) 

(0.9 x 0.3 + 0.1 x 1.2)0.0026 

/fc-snii 



I JFc-SNII \ " 
V 0.01 ) 



(3) 



For the fiducial values, the discrepancy is by a factor of 
about 6, as has been found by previous studies (e.g., Tozzi 
& Norman 2001). From Figure 1, we see that none of the 
conventional IMFs can solve the iron problem, which ranges 
from a factor 3.6 [for a Gould, Bahcall, & Flynn (1997) IMF] 
to a factor 6.6 [for a Scalo (1998) IMF] in the excess of the 
observed iron mass. From the figure, we can also see to what 
degree the IMF must be made "top heavy" in order to solve 
the problem with core-collapse SNe only; mi om needs to be 
greater than 2M Q , or a > —1.9, or some other suitable 
combination of mi ow > O.IMq and a > —2.35. 

The discrepancy can be lowered somewhat by assuming 
a larger stellar mass fraction, /„. However, only a completely 
unrealistic value of /* « 1, would lower the discrepancy to a 
non-crisis level. Stated differently, the iron mass within the 
cluster galaxies is the amount expected from core-collapse 
SNe and a normal IMF, while all the iron in the ICM is an 
excess over this expectation. This, too, has been found in 
detailed modeling (e.g. Brighenti & Mathews 1998). 

If the dominant contribution to the observed iron mass 
in clusters is from SNe la, it is straightforward to predict 
the rate of these events, Ri a , as a function of time t or 
redshift z. As already noted, the SN la rate vs time in a given 
environment is the convolution of the star-formation history 
in that environment with a "delay" or "transfer function", 
D(t), which is the SN la rate vs time following a brief burst 
of star formation at t = 0. In the case of clusters, all evidence 
is that star formation occurred considerably before z ~ 1. 
Let us assume, first, that star formation occurred in a brief 
burst, at cosmic time tf, corresponding to a redshift zf. 
Then simply 



R Ia (t) oc D(t-t f ). 



(4) 



Various authors have used different approaches to rep- 
resent D(t). For example, Ruiz-lapuente & Canal (1998) 
and Yungelson & Livio (2000) have attempted to derive 
physically motivated versions of D(t) for the different pro- 
genitor scenarios, including binary evolution and accretion 
physics. Such calculations are complex and, by necessity, 
include a large number of assumptions and poorly known 
parameters. Nevertheless, the resulting D(t) functions have 
a number of generic features: a delay until the progenitor 
population forms; a fast rise to maximum; and a power-law 
or exponential decay. These general forms suggest an alter- 
native, more phenomenological, parametrization of D(t), as 
has been adopted by Sadat et al. (1998), MDP, and Dahlen 
& Fransson (1999). 

We follow the delay function parameterization given 
by MDP. It is assumed that the progenitors of SNe la are 
white dwarfs, and therefore the overall time delay includes 
the mass-dependent lifetime of the progenitor as a main- 
sequence star, Aijy/s- Once the progenitor has become a 
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2 4 6 8 10 

Time since starburst [Gyr] 



/ 



Rla(t) dt : 



Mi 



Fc — observed 



(1- 



M Fe - 



A/- 



Fc — observed 



(7) 



The mean iron yield of a single SN la, rrtFc-ia, is gen- 
erally agreed to be about 0.7 ± 0.1M©. This emerges from 
modeling of the bolometric light curves of SNe la (e.g., Con- 
tardo, Leibundgut, & Vacca 2000, and references therein), 
as well as from SN la model calculations (e.g., Thielemann, 
Nomoto, & Yokoi 1986). 

If star formation is not in an instantaneous burst, but 
rather begins at tf and lasts some period of time, the oc- 
currence of some SNe la will be delayed. The SN la rate at 
later times, given the same total iron mass, will be neces- 
sarily greater than under the brief burst assumption. The 
assumption of a single, brief, star-formation burst at tf will 
therefore lead to a lower limit for the predicted SN la rate 
at later times. 

Rewriting Eq. 7 with the fiducial values above and in 
Eq. 3, and normalizing the SN rate by the present-day B- 
band stellar luminosity of the cluster, we obtain 



Rla(t) dt 



Figure 2. Examples of model SN la delay functions, following the 
prescription of MDP. All curves are normalized to have the same 
area, i.e., to give a time- integrated total number of SNe la that 
will produce the observed iron-mass-to-light ratio of rich clusters. 



white dwarf, it has a probability oc exp(— — ) to explode as 
a SN la, where At is the time since the star left the main 
sequence. Following MDP, then 



D(t) oc 



/ 



, t - Mms s dN 

cxp( )- — dm. 

t dm 



(5) 



For consistency, a Salpeter (1955) IMF is assumed. The min- 
imum and maximum initial masses that will lead to the for- 
mation of a WD that explodes as a SN la are 



max[3M0, ( 



t-t' 
10 Gyr^ 



M(7 



8M , 



and 
At A 



10 Gyr 



( — r 



(6) 



Figure 2 shows several examples of the MDP delay function. 
After an instantaneous starburst* at t = 0, the function is 
zero for 55 Myr (until the first white dwarfs form), then rises 
approximately as t 5 , and reaches a peak at t — 0.64 Gyr. 
It then declines exponentially with a timescale r. 

In a given cluster, the integral on R(t) over time gives 
the total number of SNe la that have exploded in that clus- 
ter. If most of the cluster iron is from SNe la, this total 
number is just the observed iron mass (minus the small frac- 
tion expected from core-collapse SNe arising from a normal 
IMF, Eq. 3), divided by the iron mass yield per SN la. The 
normalization of R(t) is therefore set by 

* See MDP, Eq. 4, for the treatment of non-instantaneous star 
formation. 



0.042. 



SN 



^ A, / 



200 



/bar \ {fgas ^Fe — gas ~\~ /* ^Fe*) 



( /bar \ 
VO.17^ 



(0.9 x 0.3 + 0.1 x 1.2)0.0026 



7n Fc -Ia 

0.7M B 



• M Fe _ obB , 



5/6 



(8) 



Here we have used a value of M/ Lb = (200 ±50) j^- for the 
typical total mass-to-light ratio measured in rich clusters, 
assuming H = 70 km s _1 Mpc -1 . The central value and 
the error is obtained by taking the union of several recent 
determinations and their quoted uncertainties - Carlberg et 
al. (1996), Girardi & Giuricin (2000), Bahcall & Comerford 
(2002), and Girardi et al. (2002). Galaxy kinematics, model- 
ing of the X-ray emission, and strong and weak gravitational 
lensing give generally consistent results for this parameter. 
The fraction of the cluster mass that is in baryons, /bar, is 
10-30% (Ettori & Fabian 1999; Mohr, Mathiesen & Evrard 
1999; Allen et al. 2002; Arnaud et al. 2002), and is thought 
to be representative of the universal baryonic mass fraction. 
We therefore adopt the cosmic value, 0.17, as recently mea- 
sured with WMAP (Spergel et al. 2003). 

The time-integrated number of SNe la per unit stellar 
luminosity obtained from Eq. 8 for the fiducial values, can 
also be expressed as 

0.042-^- = 42 SNu Gyr. (9) 

Lbq 

In other words, to produce the iron mass seen in clusters 
with SNe la, there must have been in the past one SN la for 
every 231/b© of present-day stellar luminosity. Equivalently, 
the mean SN la rate over a ~ 10 Gyr cluster age must have 
been 4.3 SNu. Since present day rates (e.g., in Virgo cluster 
ellipticals; Cappellaro et al. 1999) are much lower than this, 
the rate must have been much higher in the past. 

Thus, the assumption of a brief star-formation burst at 
some time tf in the past, followed by a SN la rate D(t — tf) 
with characteristic time r, determines the form of Ri a (t). 
The observed iron mass determines the normalization of 
Ria(t). Extended or multiple starbursts after tf can only 
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Figure 3. Predicted SN la rates vs redshift, if most of the iron 
mass in clusters is produced by type-la SNc following a brief burst 
of star formation at redshift Zf = 2 (solid curves) and Zf = 3 
(dashed curves). The different curves are for SN la transfer func- 
tions with mean delay times, r, as marked. Cluster SN la rate 
measurements are by Reiss (2000) and Gal- Yam ct al. (2002). The 
latter are shown with 95%-confidence vertical error bars. The hor- 
izontal error bars give the visibility-time- weighted redshift ranges 
of the cluster samples. 

lead to a higher Ri a {t) (except at times very soon after 
tf). The only free parameters are therefore tf and r, which 
can be constrained by comparison to direct measurements 
of Ria{t). Time and redshift are related, for our chosen cos- 
mology, by 

At = Hq 1 / (l + zy^Qmii + z) 3 + n A r°- 5 dz. (10) 



3 COMPARISON TO OBSERVATIONS AND 
CONCLUSIONS 

Figure 3 shows the expected Ri a vs. redshift based on Eqns. 
5 and 8, for several choices of r and cluster stellar formation 
redshift Zf, and compares these predictions to the observa- 
tions. Existing measurements of cluster SN la rates are by 
Reiss (2000) for 0.04 < z < 0.08, and by Gal- Yam et al. 
(2002) for z = 0.25lo.J?, and z = 0.90t£o£- The error bars 
show the 95% confidence intervals. Figure 4 is the same, but 
zooms in on the region z < 1.4, where the data exist. 

Figure 4 shows that the z ~ 1 SN la rate measurement 
is inconsistent with several of the plotted models. However, 
in the comparison, we must keep in mind that there are un- 
certainties in the parameters entering Eq. 8, which sets the 
normalization of the curves. The main uncertainty, ±25%, 
is in the mean M/Lb of clusters. Accounting also for the 
(smaller) uncertainties in the other parameters, we consider 



Figure 4. Same as Fig. 3, but zooming in on z < 1.4. The Zf = 2 
models (solid curves) with r > 2 Gyr are clearly ruled out by the 
z ~ 1 SN-rate measurement, even after accounting for a 30% 
uncertainty in the nomalization of the models. The Zf = 3 model 
(dashed curves) with t = 5 Gyr predicts unacceptably high rates 
at low 2. 

the measured 95% upper limit on Ri a at z = 0.9 to be in 
conflict with the model with Zf — 2, r = 2 Gyr, and the 
model with Zf = 2, r = 3 Gyr. The upper limit on Ri a (t) 
is at ~ 60% of the predicted values for these models. The 
model with Zf =3, r = 5 Gyr is marginally consistent with 
the z = 0.9 data point, given the uncertainty in the predic- 
tion. However, this model gives an unacceptably high rate 
at low redshift, and therefore can also be rejected. The low 
observed SN la rate at z ~ — 1 means that, if cluster iron 
was produced by SNe la, those SNe must have occurred at 
earlier times, times that have yet to be probed by obser- 
vations. To push the SNe la to such early times, one must 
invoke early star formation and a short SN la delay time. 

A recent attempt by van Dokkum & Franx (2001) to 
deduce the epoch of star formation in cluster ellipticals by 
spectral synthesis modeling of the stellar populations, and 
accounting for biases in the selection of elliptical galaxies, 
finds stellar formation redshifts of Zf = 2+ ' 2 . If we adopt 
a particular formation redshift, say Zf = 2, the resulting 
upper limit of r < 2 Gyr places clear constraints on SN la 
progenitor models. Two of the models by Yungelson & Livio 
(2000) predict longer delays and are therefore ruled out. In 
the double degenerate model, two WDs merge, and the re- 
sulting object, having a mass larger than the Chandrasekhar 
limit, is subject to a runaway thermonuclear explosion trig- 
gered at its core. Yungelson & Livio find that the delay 
function of such systems has an exponential cutoff at ~ 11 
Gyr, at odds with the above upper limit on r. A similar 
discrepancy exists for a model where a WD accrets He-rich 
material from a He-star companion, leading to helium igni- 
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tion on the surface of the WD and an edge-lit detonation 
of the star. This model has a delay function that cuts off at 
~ 5 Gyr. 

One must recap here that all these conclusions hold 
under the assumption that the stars in galaxy clusters were 
formed with a standard IMF, and therefore most of the iron 
is from SNe la. If the stars were formed with a sufficiently 
top-heavy IMF to produce the observed iron, few SNe la are 
expected at any redshift, and the SN la rate measurements 
place no constraints on cluster star-formation epoch or on 
SN la time delay. 

Since our conclusions are based on the comparison of a 
low observed cluster SN la rate at z ~ 1 to the high rate 
predicted by some models, it is sensible to re-examine the 
reliability of the observation. One possibility to consider is 
that the rate measured by Gal- Yam et al. (2002), based on 
deep HST cluster images, was low because the SN detec- 
tion efficiency was overestimated. Some faint cluster SNe 
that were missed would then be incorrectly accounted for in 
the rate calculations. This is highly unlikely, given that the 
actually detected cluster SNe were relatively bright, with 
I < 24 mag. Much fainter SNe, down to 28 mag, were found 
by the survey, demonstrating its high sensitivity, but these 
SNe were background and foreground events, rather than 
cluster events. A second possibility is that the rate measured 
was skewed low by considering only the visibility time of 
normal SNe la. Ignoring a significant population of sublumi- 
nous SNe la, which have short visibility times, will lower the 
derived rate. However, studies of local subluminous SNe la 
show that they have an extremely low iron yield, ~ O.OO7M 
(e.g., Contardo et al. 2000). Thus, unless they are very com- 
mon at z ~ 1, such SNe are irrelevant for iron production in 
clusters. 

An interesting conclusion arises if we examine our re- 
sults jointly with the results we report in a companion paper 
(Gal- Yam & Maoz 2003). There, we find that some versions 
of cosmic SFH, combined with particular SN la delay times, 
are incompatible with the observed redshift distribution of 
SNe la found by Perlmutter et al. (1999). For example, if 
the cosmic SFR rises between z = and z ~ 1 as sharply as 
implied by Lilly et al. (1996) and by Hippelein et al. (2003), 
and at z > 1 as sharply as found by Lanzetta et al. (2002), 
then the observed SN la redshift distribution sets a 95%- 
confidence lower limit on the SN la delay time, of r > 3 Gyr. 
Thus, if one had independent reasons to adopt this SFH sce- 
nario, and a cluster star- formation redshift Zf < 2 (the lat- 
ter implying an upper limit of r < 2 Gyr at 95% confidence 
from the measured SN la rate), then all values of r would be 
excluded. One would then be forced to the conclusion that 
SNe la cannot be the source of iron in clusters, leaving only 
the top-heavy IMF option. 

To summarize, we have investigated the source of the 
total iron mass in rich galaxy clusters. Using updated val- 
ues for the various observational parameters, we have red- 
erived and quantified the excess of iron over the expectation 
from core-collapse SNe, provided the stars in cluster galax- 
ies formed with a standard IMF. Assuming the source of the 
iron is from SNe la, we then showed that the SN la rate vs. 
redshift can be predicted quite robustly, given two parame- 
ters - the star formation redshift in clusters, and the mean 
SN la delay time. We set constraints on these two parame- 
ters by comparing the predicted rates to the measurements 



of the SN la rate at z ~ 1 by Gal- Yam et al. (2002). The 
low observed rates at z ~ 1 force the iron-producing SNe la 
to have occurred at higher redshifts. This implies an early 
epoch of star formation and a short delay time. Specifically, 
we showed that models with a mean SN la delay time of 
r < 2 Gyr (< 5 Gyr) and Zf < 2 (zf < 3) are ruled out at 
high confidence. Thus, if other avenues of inquiry show that 
SNe la explode via a mechanism that leads to a long delay 
time that is ruled out by our study, it will mean that core- 
collapse SNe from a top-heavy IMF must have formed the 
bulk of the observed mass of iron in clusters. The same con- 
clusion will apply if other studies confirm a cosmic SFH that 
rises sharply with redshift. As shown in Gal- Yam & Maoz 
(2003), short (< 4 Gyr) delay times are incompatible with 
such a SFH and the redshift distribution of field SNe la. 

Finally, we note that our results have relied on the up- 
per limits on the z ~ 1 cluster SN la rate, calculated by 
Gal- Yam et al. (2002) based on just several SNe. Improved 
measurements of cluster SN rates at low, intermediate, and 
high redshifts can tighten the constraints significantly, and 
can potentially reveal the iron source directly. 
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